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Abstract 



Event-by-event multiplicity fluctuations in central C+C, S+S, In+In, and Pb+Pb as well as p+p 



collisions at bombarding energies from 10 to 160 AGeV are studied within the HSD and UrQMD 
microscopic transport approaches. Our investigation is directly related to the future experimental 
program of the NA61 Collaboration at the SPS for a search of the QCD critical point. The dependence 
on energy and atomic mass number of the scaled variances for negative, positive, and all charged 

(N 

hadrons is presented and compared to the results of the model of independent sources. Furthermore, 

""xl" ! 

the nucleus-nucleus results from the transport calculations are compared to inelastic proton-proton 
collisions for reference. We find a dominant role of the participant number fluctuations in nucleus- 
nucleus reactions at finite impact parameter b. In order to reduce the influence of the participant 

; 

numbers fluctuations on the charged particle multiplicity fluctuations only the most central events 



have to be selected. Accordingly, the samples of the 1% most central nucleus-nucleus collisions with 
the largest numbers of the projectile participants are studied. The results are compared with those 
for collisions at zero impact parameter. A strong influence of the centrality selection criteria on the 
multiplicity fluctuations is pointed out. Our findings are essential for an optimal choice of colliding 
nuclei and bombarding energies for the experimental search of the QCD critical point. 

PACS numbers: 24.10.Lx, 24.60.Ky, 25.75.-q 
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I. INTRODUCTION 



The event-by-event fluctuations in high energy nucleus-nucleus (A+A) collisions (see e.g., 
the reviews are expected to be closely related to the transitions between different 

phases of QCD matter. By measuring the fluctuations one should observe anomalies from 
the onset of deconfinement Q] and dynamical instabilities when the expanding system goes 
through the 1-st order transition line between the quark-gluon plasma (QGP) and the hadron 
gas 5[. Furthermore, the QCD critical point may be signaled by a characteristic pattern in 
enhanced fluctuations. A+A collisions in the SPS energy region are expected to be a suitable 
tool for a search of critical point signatures [6j. Only recently first measurements of particle 
multiplicity fluctuations and transverse momentum fluctuations jij in A+A collisions have 
been performed. A theoretical analysis of multiplicity fluctuations for the hadron-resonance 
gas - in different statistical ensembles - has been performed in Ref. 9]. Independently, the 
multiplicity fluctuations in A+A collisions has been studied within a microscopic transport 



approach 
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111 ] . We recall that fluctuations traditionally are quantified by the ratio of the 



variance of the multiplicity distribution to its mean value, the scaled variance. Previous works 



on this subject have to be quoted. The calculations of the statistical models [13| and transport 
approaches HSD 141 and UrQMD v. 1.3 [jjj] have been compared with the corresponding 
preliminary results 16J of the NA49 Collaboration in central Pb+Pb collisions at SPS energies. 
At RHIC energies a HSD analysis of the preliminary data [uj of the PHENIX Collaboration 
in Au+Au collisions at y/s^N = 200 GeV has been presented in Ref. [181 ]. 

An ambitious experimental program for the search of the QCD critical point has been started 
by the NA61 Collaboration at the SPS [fj]]. The program includes a variation in the atomic 
mass number A of the colliding nuclei as well as an energy scan. This allows to scan the phase 
diagram in the plane of temperature T and baryon chemical potential \xb near the critical point 
as argued in Ref. [19|. One expects to 'locate' the position of the critical point by studying 
its 'fluctuation signals'. High statistics multiplicity fluctuation data will be taken for p+p, 
C+C, S+S, In+In, and Pb+Pb collisions at bombarding energies of £7^=10, 20, 30, 40, 80, 
and 158 AGeV. 

The aim of the present paper is to study the energy and system size dependence of event- 
by-event multiplicity fluctuations within the microscopic transport approaches Hadron-String- 
Dynamics (HSD, v. 2.5) (20}] and Ultra-Relativistic-Quantum-Molecular-Dynamics (UrQMD, 
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v 1.3) [2lJ. These models provide a rather reliable description (see, e.g., Refs. |2jJ, |22| ) for 
the inclusive spectra of charged hadrons in A+A collisions from SIS to RHIC energies. In our 
study we will consider C+C, S+S, In+In, and Pb+Pb collisions at the bombarding energies 
of 10, 20, 30, 40, 80, 158 AGeV. For a comparison and reference we also present the 
results of multiplicity fluctuations in p+p collisions at the same energies. Our study thus is in 
full correspondence to the experimental program of the NA61 Collaboration [lj]]. In order to 
estimate systematical errors from the theoretical side we employ the two independent transport 
models (HSD and UrQMD) as in Ref. |22|. 

The QCD critical point is expected to be experimentally seen as a non-monotonic dependence 
of the multiplicity fluctuations, i.e. a specific combination of atomic mass number A and 
bombarding energy Ei a b could move the chemical freeze-out of the system close to the critical 
point and show a 'spike' in the multiplicity fluctuations. Since HSD and UrQMD do not 
include explicitly a phase transition from a hadronic to a partonic phase, we can not make a 
clear suggestion for the location of the critical point - it is beyond the scope of our hadron- 
string models. However, our study might be helpful in the interpretation of the upcoming 
experimental data since it will allow to subtract simple dynamical and geometrical effects from 
the expected QGP signal. The deviations of the future experimental data from the HSD and 
UrQMD predictions may be considered as an indication for the critical point signals. 

Theoretical estimates give about 10% increase of the multiplicity fluctuations due to the 
critical point 6|. It is large enough to be observed experimentally within the statistics of 
NA61 [19|. To achieve this goal, it is necessary to have a control on other possible sources 
of fluctuations. One of such sources is the fluctuation of the number of nucleon participants. 
It has been found in Ref. |l0|] that these fluctuations give a dominant contribution to hadron 
multiplicity fluctuations in A+A collisions. On the other hand one can sup pre ss the participant 
number fluctuations by selecting most central A+A collisions (see Ref. 10, 2^] for details). 
That's why the NA61 Collaboration plans to measure central collisions of light and intermediate 
ions instead of peripheral Pb+Pb collisions. It is important to stress, that the conditions for 
the centrality selection in the measurement of fluctuations are much more stringent than those 
for mean multiplicity measurements. This issue will be discussed in detail in our paper. 

Our paper is organized as follows. In Section II the HSD and UrQMD models are compared 
with the data for the charged hadron multiplicity, mean values and fluctuations in p+p collisions 
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in the SPS energy range 10 — 158 AGeV. In Section III the participant number fluctuations 
in fixed target experiments are discussed and estimated within HSD and UrQMD. In Section 
IV we present the HSD and UrQMD results for the charged hadron multiplicity fluctuations in 
A + A collisions with zero impact parameter, 6 = 0. The participant number fluctuations in 
A + A collisions at b = are considered in Section V. In addition, the transport model results 
for the charged hadron multiplicity fluctuations are compared within the model of independent 
sources. In Section VI the centrality selection - by fixing the number of projectile participants 
- is considered and compared to the case of b = 0. Section VI presents also the HSD and 
UrQMD results for the charged hadron multiplicity fluctuations in A + A collisions for the 
1% most central collisions corresponding to the largest number of projectile participants. A 
summary and conclusion will close the paper in Section VII. 



II. MULTIPLICITY FLUCTUATIONS IN PROTON-PROTON COLLISIONS 

For a quantitative measure of the particle number fluctuations it is convenient to use the 
scaled variances, 

w - = — m — • (1) 

where (••■) denotes event-by-event averaging and the index i means "+", and "ch", i.e 
negative, positive, and all charged final state hadrons. 

The energy dependence of the measured charged multiplicity and fluctuations for p+p col- 
lisions can be parametrized by the functions [1]: 



(N ch ) - -4.2 + 4.69 (^)°' 31 , * ch = 0.35 (( ^~ 1)2 , (2) 

where ^Snn is the center-of-mass energy. 

Figure [1] shows the HSD and UrQMD results for inelastic p+p collisions in comparison to 
the experimental data taken from Ref. jl|. As seen from Fig. [TJ both models give a good 
reproduction of the p+p data for (N c h) , but slightly (over) underestimate u c h at high collision 
energies. The differences between the HSD and UrQMD model results for u c h can be attributed 
to different realizations of the string fragmentation model, in particular, differences in the 
fragmentation functions and the fragmentation scheme, i.e. fragmentation via heavy baryonic 



and mesonic resonances in UrQMD and direct light hadron production by string fragmentation 
in HSD. 




E (GeV) E (GeV) 

FIG. 1: (Color online) The average multiplicity (left) and scaled variance (right) of charged hadrons 
in p+p inelastic collisions. The open circles and squares (connected by solid lines) show the results of 
HSD and UrQMD, respectively, whereas the full circles present the experimental data from Ref. 



For negative and positive charged hadrons the average multiplicities and scaled variances 
in p+p collisions can be presented in terms of the corresponding quantities for all charged 
particles, 

<JV±> = \ ((Nch) ± 2) , ^=l^^giL_. (3) 

III. PARTICIPANT NUMBER FLUCTUATIONS 

In each A + A collision only a fraction of all 2A nucleons, i.e. the participants, interact. 
Participants from the projectile and target nuclei are denoted as N p r ° 3 and Np Tg , respectively. 
The nucleons, which do not interact, are denoted as projectile and target spectators, Ng ro3 = 
A — Np OJ and N^' 9 = A — N 1 ^ 9 . The spectators are selected according to the criteria \y — 
ybeam(target)\ < 0.32 and excluded from the multiplicity fluctuation analysis. We recall that 
fluctuations in high energy A + A collisions are dominated by a geometrical variation of the 
impact parameter b. However, even for fixed impact parameter b the number of participants, 
N P = Np ro:) + Np r9 , fluctuates from event to event. This is due to fluctuations in the initial 
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states of the colliding nuclei and the probabilistic character of the interaction process. These 
fluctuations of Np form usually a large and uninteresting background. 

To minimize the event-by-event fluctuations of the number of nucleon participants in measur- 
ing the multiplicity fluctuations the NA49 Collaboration has been trying to fix N p roj in Pb+Pb 
collisions. Samples of collisions with a fixed number of projectile spectators, Ng r ° 3 = const 
(and thus a fixed number of projectile participants, N p r ° 3 ), have been selected. This selection 
is possible in fixed target experiments at the SPS, where Ng TOJ is measured by a Zero Degree 
Veto Calorimeter covering the projectile fragmentation domain. A similar centrality selection 
is expected to be implemented in the future NA61 experiment. However, even in samples 
with N p roj = const the number of target participants will fluctuate considerably. Hence, an 
asymmetry between projectile and target participants is introduced, i.e. N p r ° 3 is constant by 
constraint, whereas Np 0,1 ' 9 fluctuates independently (the consequences of this asymmetry have 
been discussed in Ref. |23|). 

In each sample with N p r ° 3 = const the number of target participants fluctuates around 
its mean value with the scaled variance u/p r9 . The mean value equals to (Np ar9 ) = N p r ° 3 , if 
N p r ° 3 is not too close to its limiting values, N p r ° 3 = 1 and N p r ° 3 =A. The scaled variance of 
target participants U3 l p 9 as a function of fixed number of projectile participants N p r ° 3 has been 
obtained from HSD and UrQMD in Ref. (loj for Pb+Pb collisions at 158 AGeV. 

Fig. [2] (left) presents the HSD scaled variances ui p r9 for C+C, S+S, In+In, and Pb+Pb 
collisions at 158 AGeV as a function of N p ° 3 . The fluctuations of N^ 9 are quite strong for 
peripheral reactions (small N p roj ) and negligible for the most central collisions (large N p r ° 3 ). 
A vanishing of u t p r9 = at N p r ° 3 =A does not, however, show up in C+C collisions. Even for 
N p p ° 3 =A=12 in C+C collisions the number of participants from the target still fluctuates and 
the scaled variance amounts to <jj l P 9 = 0.25. Fig. [2] (right) shows ufp 1 " 9 for light nuclei, S+S, 
Ne+Ne, O+O, and C+C. Even for the maximal values of N p ° 3 =A the fluctuations uf^ 9 do 
not vanish and increase with decreasing atomic mass number A. 

The temperature T and baryon chemical potential fip at the hadron chemical freeze-out 
demonstrate the dependence on both the collision energy and system size 12J. Thus, changing 
the number of participating nucleons one may scan the T — fip plane. Some combination of 
Np and Ei a b might move the chemical freeze-out point close the QCD critical point. One could 
then expect an increase of multiplicity fluctuations in comparison to their 'background values'. 
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FIG. 2: (Color online) Left: The scaled variance uip rg for the fluctuations of the number of target 
participants, Np rg . The HSD simulations of u}p T9 as a function of N p ro:i are shown for different 
colliding nuclei, In+In, S+S, and C+C at £ /a6 =158 AGeV. The HSD and UrQMD results for Pb+Pb 
collisions are taken from Ref. [lo| . Right: The scaled variance Wp r9 for light nuclei, S+S, Ne+Ne, 
O+O, and C+C at Ei a b=158 AGeV as a function of N^ 03 . Fluctuations of target participants for 
jyproj _^ _ s } lown by the open symbols - are different from zero and increase with decreasing atomic 
mass number A. 

Why does one need central collisions of light and intermediate ions instead of studying 
peripheral Pb+Pb collisions for a search of the critical point? Fig. [2] explains this issue. At fixed 
N p r ° 3 the average total number of participants, Np = N p r ° 3 + N^ 9 , is equal to (Np) = 2N p roj , 
and, thus, it fluctuates as uj p = 0.5LUp T9 . Then, for example, the value of Np™ 3 = 30 corresponds 
to almost zero participant number fluctuations, Up = 0, in S+S collisions while up becomes 
large and is close to 1 and 1.5 for In+In and Pb+Pb, respectively. Even if N p roj is fixed exactly, 
the sample of the peripheral collision events in the heavy-ion case contains large fluctuations 
of the participant number: this would 'mask' the critical point signals. As also seen in Fig. [2] 
(right), the picture becomes actually more complicated if the atomic mass number A is too 
small. In this case, the number of participants from a target starts to fluctuate significantly 
even for the largest and fixed value of N p ro:! =A. 



7 



IV. MULTIPLICITY FLUCTUATIONS AT ZERO IMPACT PARAMETER 



A. Centrality Selection in A+A Collisions by Impact Parameter 

The importance of a selection of the most central collisions for studies of hadron multiplicity 
fluctuations has been stressed in our previous papers [lQ, H, Q, Q. Due to its convenience 
in theoretical studies (e.g., in hydrodynamical models) one commonly uses the condition on 
impact parameter 6, for the selection of the 'most central' collisions in model calculations. 
However, the number of participant even at 6 = is not strictly fixed, and fluctuates according 



to some distributions (cf. Fig. 14 (right.) from Ref. [ill]). It should be stressed again that 
the conditions b < b max can not be fixed experimentally since the impact parameter itself can 
not be measured in a straightforward way. Actually, in experiments one accounts for the 1%, 
2% etc. most central events selected by the measurement of spectators in the Veto calorimeter, 
which corresponds to the event class with the largest Np" 3 . As we will demonstrate below the 
multiplicity fluctuations are very sensitive to the centrality selection criteria. In particular, the 
transport model results for 6 = and for 1% events with the largest N p roj are rather different 
(see below). 

Let's start with the 6 = centrality selection criterium. We recall that the charged multi- 

nn 

plicity fluctuations are closely related to the fluctuations of the number of participants [10, H4j . 
Therefore, it is useful to estimate the average number of participants, (Np), and the scaled 
variances of its fluctuations, up, in A + A collision events which satisfy the 6 = condition. 
The left panel in Fig. [3] shows the ratio, (Np)/2A, in A + A collisions with 6 = for differ- 
ent nuclei at collision energies Ei^ = 10 and 158 AGeV. Both transport models (HSD and 
UrQMD) show a monotonous increase of (Np) /2A with collision energy for all nuclei in the 
energy range 10 + 158 AGeV (Fig. [3l left). Correspondingly, the fluctuations of the number of 
participants Up for all nuclei become smaller with increasing collision energy (Fig. [31 right.). 
As seen from Fig. [3] (left) about 90% of nucleons are participants for Pb+Pb collisions with 
6 = 0. This number becomes essentially smaller, about 60-70%, for C+C collisions. One can 
therefore expect that participant number fluctuations at 6 = are small for heavy nuclei but 
strongly increase for light systems. This is demonstrated in Fig. [3] (right): up is about 0.1-i-0.2 
in Pb+Pb and In+In but becomes much larger, 0.5^0.7, in C+C collisions. 

One can conclude that the condition 6 = corresponds to 'most central' A + A collisions 
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FIG. 3: (Color online) Left: Mean (Np), divided by the maximum number of participants 2 A in events 
with 6 = for different nuclei at collision energies £^=10 and 158 AGeV. Right: The scaled variance 
u>p in events with 6 = for different nuclei at collision energies Ei a b=W and 158 AGeV. 

only for nuclei with large atomic mass number (In and Pb). In this case the average number 
of participants is close to its maximum value and its fluctuations are rather small. However, in 
the studies of event-by-event multiplicity fluctuations in the collisions of light nuclei (C and S) 
the criterium b = is far from selecting the 'most central' A + A collisions. 



B. HSD and UrQMD Results for the Multiplicity Fluctuations for 6=0 

Results of HSD and UrQMD transport model calculations for the scaled variance of negative, 
positive, u> + , and all charged, u> c h, hadrons are shown in Figs. H] and [5] at different collision 
energies, E lab = 10, 20, 30, 40, 80, 158 AGeV, and for different colliding nuclei, C+C, S+S, 
In+In, Pb+Pb. The transport model results correspond to collision events for zero impact 
parameter, 6 = 0. To make the picture more complete, the transport model results for inelastic 
p+p collisions are shown too, for reference. Note that in our presentation - throughout the 
paper - the proton spectators are not accounted for in the calculation of N + and N c h- Thus, 
proton spectators do not contribute to u + and uj c h. 

One sees a monotonic dependence of the multiplicity fluctuations on both and A: the 
scaled variances u + , and u ch increase with E\ ah and decrease with A. The results for 
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FIG. 4: (Color online) The results of HSD (left) and UrQMD (right) simulations for uj- (top panel), 
uj + (middle panel), and uj c h (lower panel) in p+p and central C+C, S+S, In+In, Pb+Pb collisions at 
Eiab = 10, 20, 30, 40, 80, 158 AGeV. The condition b = is used here as a criterium for centrality 
selection. There are no cuts in acceptance. 
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FIG. 5: (Color online) The same as in Fig. HI but only hadrons with positive cm. rapidities, y > 
(projectile hemisphere), are accepted. 
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p+p collisions are different from those for light ions. We note that within HSD and UrQMD 
a detailed comparison of the multiplicity fluctuations in nucleon-nucleon inelastic collisions 
and 6 = heavy-ion collisions (Pb+Pb and Au+ Au) , including the energy dependence up to 



y/SNN = 200 GeV, has been presented in Refs. [14j, Il5l |. 



Fig. H] corresponds to the full An acceptance, i.e. all particles are accepted without any cuts 
in phase space. In actual experiments the detectors accept charged hadrons in limited regions 
of momentum space. Fig. [5] shows the HSD and UrQMD results for multiplicity fluctuations in 
the projectile hemisphere (i.e. positive rapidities, y > in the cm. frame). This corresponds to 
the maximal possible acceptance, up to 50% of all charged particles, by the optimized detectors 



of the NA61 Collaboration 



191 ] . One observes from Fig. [5] that the energy and system size 



dependencies of the multiplicity fluctuations in the projectile hemisphere (y > 0) become less 
pronounced than in full An acceptance. Note also that the centrality selection criterium 6 = 
keeps the symmetry between the projectile and target hemispheres. Thus, the results for a 
y < acceptance are identical to those for y > presented in Fig. [5j 



C. Comparison to the Independent Source Model 

The multiplicity fluctuations in elementary nucleon-nucleon collisions and fluctuations of the 
number of nucleon participants are presented in the right panels of Figs. [1] and [3l respectively. 
Their combination explains the main features of hadron multiplicity fluctuations in A + A 
collisions shown in Figs. H] and [5J in particular, the dependence on collision energy and atomic 
mass number. They also are responsible for the larger values of ioi in the UrQMD simulations 
in comparison to those from HSD. To illustrate this let us consider the model of independent 
sources (ISM). 

The multiplicity fluctuations in A + A collisions can be then written according to the ISM 



as (see e.g., Refs. [l 







14 Ha, 123), 

Ui = U* + Hi Up , (4) 



where u* denotes the fluctuations of the hadron multiplicity from one source and the term rij up 
gives additional fluctuations due to the fluctuations of the number of sources. One usually 
assumes that the number of sources is proportional to the number of nucleon participants. 
The value of Hi in Eq. (J4j) then is the average number of z'th particles per participant, rij = 

12 



(Ni) / (Np), and ujp equals the scaled variance for the number of nucleon participants. Nucleon- 
nucleon collisions, which are the weighted combinations of p+p, p+n, and n+n reactions, define 
the fluctuations uj* from a single source (see details in Refs. jlO, Q). 




j i i i i i i i i i i i i i i i l i ... i i 



20 40 60 80 100 120 140 160 40 80 120 160 200 

E (GeV/nucleon) A, mass number 

FIG. 6: (Color online) The left panel illustrates the energy dependence of lo- in S+S collisions at 
b = in the full 4ir acceptance, the right panel - the uo- dependence on atomic mass number at 
Ei a b = 80 AGeV. The HSD results are shown by the squares while the circles correspond to Eq. (f4]) of 
the model of independent sources. The stars show the first term, uj*_, in the right, of Eq. (T4|) - scaled 
variance for negative hadrons in nucleon- nucleon collisions, where the values of uj*_, nj, and ojp are 
calculated within HSD. 

In Fig. O the HSD results for uj^ in A + A collisions at b = are compared to the ISM 
- Eq. P|. One concludes that the transport model results for the multiplicity fluctuations 
are in qualitative agreement with Eq. (j3J) of the independent source model. Both m and uj* 
increase strongly with collision energy as seen from Fig. UJ This explains, due to Eq. (j3J), the 
monotonous increase with energy of the scaled variances uji in A + A collisions at b = seen in 
Figs. H] and El Note that up at b = decreases with collision energy as shown in Fig. [21 right 
This, however, does not compensate a strong increase of both Uj and uj*. The atomic mass 
number dependence of the scaled variances uji in A + A collisions with 6 = follows from the 
A-dependence of ujp. Fig. [3] (right) demonstrates a strong increase of ujp for light nuclei. This, 
due to Eq. is transformed to the corresponding behavior of uji seen in Figs. H] and [51 
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V. MULTIPLICITY FLUCTUATIONS IN 1% MOST CENTRAL COLLISIONS 



We consider now the centrality selection procedure by fixing the number of projectile par- 
ticipants N p ro:> . This corresponds to the real situation of A + A collisions in fixed target 
experiments. As a first step we simulate in HSD and UrQMD the minimal bias events - which 
correspond to an all impact parameter sample - and calculate the event distribution over the 
number of participants N part . Then, we select 1% most central collisions which correspond to 
the largest values of N^ OJ . In such a sample of A+A collisions events with largest N'p 03 from 
different impact parameters can contribute. After that we calculate the values of up in these 
samples. Note that even for a fixed number of N p ro:) the number of target participants N t p Tg 
fluctuates. Thus, the total number of participants, Np = N l p 79 + N p ro3 , fluctuates too. In our 
1% sample, both Np r9 and N^° 3 fluctuate. Besides, there are correlations between Np r9 and 

n p ;° 3 . 

Fig. [7] shows the ratio {Np)/2A and the scaled variance, up, for 1% most central collisions 
selected by the largest values of N^" 3 . These results are compared with those for the 6 = 
centrality selection. For heavy nuclei, like In and Pb, one finds no essential differences between 
these two criteria of centrality selection. However, the 1% centrality trigger defined by the 
largest values of N p ro:i looks much more rigid for light ions (S and C). In this case, the ratio 
(Np)/2A is larger, and uop is essentially smaller than for the criterion 6 = 0. As a result, the 
1% centrality trigger by the largest values of N p ro:> leads to a rather weak A-dependence of top. 

Some comments are appropriate at this point. Let us define the centrality c(N) as a percent- 
age of events with a multiplicity larger than N (this can be the number of produced hadrons, 



number of participants, etc.). It was argued in Ref. 24| that a selection of c(N) of most central 



A + A collisions is equivalent to restricting the impact parameter, b < b(N), with, 



b(N) = J^f- c(N) , (5) 

where <Ji ne i is the total inelastic A+A cross section. Thus, the centrality criterion by the 
multiplicity N is equivalent to the geometrical criterion by the impact parameter b. Moreover, 
the result (jSJ) does not depend on the specific observable N used to define the c-percentage 
of most central A+A collisions. Eq. (jSJ) should remain the same for any observable N which 
is a monotonic function of b. Therefore, the relation (0) reduces any centrality selection to 



the geometrical one. This result was obtained in Ref. 
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24] by neglecting the fluctuations of 
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FIG. 7: (Color online) The HSD (left) and UrQMD (right) results for the ratio (N P )/2A (the upper 
panel) and the scaled variance of the participant number fluctuations, top (the lower panel), for the 
1% most central collisions selected by the largest values of Np TO:! (full symbols), for different nuclei at 
collision energies £'; a b=10 and 158 AGeV. The open symbols present the results of Fig. [3] (right) for 
6 = 0. 

multiplicity N at a given value of b. This is valid if c is not too small and the colliding nuclei 
are not too light. In the sample of A+A events with 1% of largest N p ro:> , the relation (jSJ) can 
not be applied for S+S and C+C collisions. The average value of (N p ro:1 ) even for b = is 
essentially smaller than its maximal value A. To form the sample with 1% largest N"^° 3 , one 
needs at several fixed values of b to take into account the fluctuations with N p r ° 3 > (Np ro3 (b)). 
Just the fluctuations of the number of nucleon participants form the 1% sample with largest 
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FIG. 8: (Color online) The dependence of uj- on atomic mass number at £^=80 AGeV for the HSD 
(left) and UrQMD (left) simulations. The squares correspond to b = 0, and circles to 1% largest Npf° 3 . 

Np° 3 values. 

Fig. [8] shows a comparison of the A-dependence of uj- in the transport models for two 
different samples of the collision events: for 6 = and for the 1% of events with largest N^° 3 
values. One can see that the multiplicity fluctuations are rather different in these two samples. 
Moreover, these differences are in the opposite directions for heavy nuclei and for light nuclei. 
For light nuclei, u>- is essentially smaller in the 1% sample with largest N"p° 3 values, whereas 
for heavy nuclei the smaller fluctuations correspond to 6 = events. Note that in the 1% 
sample with largest N[p° 3 values the A-dependence of multiplicity fluctuations becomes much 
weaker. In this case, a strong increase of the multiplicity fluctuations for light nuclei, seen for 
6 = 0, disappears. 

For the 1% most central A + A collision events - selected by the largest values of N^° 3 - 
the HSD multiplicity fluctuations are shown in Figs. [9] and [11] and the corresponding UrQMD 
results are shown in Figs. [10] and [121 The results from both models are also presented in Tabled] 
The model uncertainties are shown as errorbars in Figs. [9] and [10] For light nuclei (S and C) 
the multiplicity fluctuations in the samples of 1% most central collisions are smaller than in the 
6 = selection and the atomic mass number dependencies become less pronounced (compare 
Figs. [9] and [10] with Fig. H]). This is because the participant number fluctuations u P have now 
essentially smaller A-dependence, as seen in Fig. 
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FIG. 9: (Color online) The HSD results for u;_ (upper panel), u> + (middle panel), and u} c h (lower 
panel) in A + A and p+p collisions for the full 4tt acceptance in 3D (left) and 2D (right.) projection. 
The 1% most central C+C, S+S, In+In, and Pb+Pb collisions are selected by choosing the largest 
values of Np roj at different collision energies Ei a b=W, 20, 30, 40, 80, 158 AGeV. The errorbars indicate 
the estimated uncertainties in the model calculations. The HSD results from inelastic p+p collisions 




FIG. 10: (Color online) The same as in Fig. but for the UrQMD. 
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TABLE I: The HSD and UrQMD scaled variances and uj c h for the 1% of most central collisions 

selected by largest values of N^° 3 . The numbers correspond to those presented in Figs. I9l [T2l 
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FIG. 11: (Color online) The same as in Fig. El but for final hadrons accepted in the projectile 
hemisphere, y > (left), and in the target hemisphere, y < (right). The HSD results in inelastic 
p+p collisions are the same as in Fig. [5j 
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FIG. 12: (Color online) The same as in Fig. HH but for the UrQMD. 
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Fig. H] shows that both HSD and UrQMD predict a monotonic dependence of the charge 
particle multiplicity with energy. So, the hadronic 'background' for the NA61 experiments is 
expected to be a smooth monotonic function of beam energy. 

Besides of differences in the realization of the string fragmentation model in HSD and 
UrQMD 1.3 mentioned above (cf. Fig. 1), additional deviations can be attributed to dif- 
ferent initializations of the nuclei in both models. Indeed, the event-by-event observables show 
a higher sensitivity to the initial nucleon density distribution than the standard single particle 



observables 



251 ] . A pilot study using UrQMD shows different u when applying different initial- 



ization shapes. Due to this effects a systematic error of 20% is attributed to uo. Such a sensitivity 
of the A-dependence of u) c h to the details of the models indicates a necessity for further studies 
of the initializations of the nuclei in transport model approaches. This becomes important for 
the theoretical interpretation of future experimental data on event-by-event fluctuations. 

Note that the model of independent sources and Eq. (jlj) work for the multiplicity fluctuations 
simulated by the transport models in full An acceptance but not for the acceptance in a specific 
rapidity region. The results for inelastic p+p collisions are identical in the projectile and target 
hemispheres. This is not the case in the sample of 1% most central A + A collisions selected 
by N p ro:> . The total number of nucleons participating in A + A collisions fluctuates. These 
fluctuations are not symmetric in forward-backward hemispheres: in the selected 1% sample 
the number of target participants N l p r9 fluctuates essentially stronger than that of N^° 3 . The 
consequences of the asymmetry in an event selection depend on the dynamics of A+A collision 
(see Ref. 23j for details). The HSD and UrQMD results in Figs. [TT1 and [121 clearly demonstrate 



larger values for all scaled variances, cj_, u + , and u c h, for y < acceptance than those for y > 
one. This is due to stronger target participant fluctuations, uf^ 9 > uj P p° 3 . 



VI. SUMMARY AND CONCLUSIONS 



In summary, the event-by-event multiplicity fluctuations in nucleus-nucleus collisions have 
been studied for different energies and system sizes within the HSD and UrQMD v. 1.3 transport 
approaches. Our present study is in full correspondence to the future experimental program of 
the NA61 Collaboration at the SPS. Thus we have considered C+C, S+S, In+In, and Pb+Pb 
nuclear collisions from Ei a f,= 10, 20, 30, 40, 80, 158 AGeV. The influence of participant 
number fluctuations on hadron multiplicity fluctuations has been emphasized and studied in 
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detail. To make these 'trivial' fluctuations smaller, one has to consider the most central colli- 
sions. Indeed, one needs to make a very rigid selection - 1% or smaller - of the 'most central' 
collision events. In addition, one wants to compare the event-by-event fluctuations in these 
'most central' collisions for heavy and for light nuclei. Under these new requirements different 
centrality selections are not equivalent to each other. As a consequence, there is no universal 
geometrical selection by the impact parameter. This is a new and serious problem for theoret- 
ical models (e.g., for hydrodynamical models) in a precision description of the event-by-event 
fluctuation data. The above statements have been illustrated by the 6 = selection criterium 
considered in our paper. For light nuclei even these 'absolutely central' geometrical collisions 
lead to rather large fluctuations of the number of participants, essentially larger than in the 1% 
most central collisions selected by the largest values of the projectile participants N^ OJ . 

We have, futhermore, used the number of projectile participants to define the centrality 
selection. This is the most promising way of the centrality selection in fixed target experiments. 
It also corresponds to the experimental plans of the NA61 collaboration. We have defined the 
1% most central collisions by selecting the largest values of the projectile participants N^° 3 . 
The multiplicity fluctuations calculated in these samples show a much weaker dependence on 
the atomic mass number A than for criterium b = 0. A monotonic energy dependence for the 
multiplicity fluctuations are obtained in both the HSD and UrQMD transport models. The 
two models demonstrate a similar qualitative behavior of the particle number fluctuations. 
However, the UrQMD 1.3 results for the scaled variances cj_, u>+, and u) c h are systematically 
larger than those obtained within HSD. This is mainly due to the corresponding inequalities 
for the scaled variances uj c h (see Fig. 1, right) for p+p collisions in these models. Our study 
has demonstrated a sensitivity of the multiplicity fluctuations to some specific details of the 
transport models. Nevertheless, the present HSD and UrQMD results for the scaled variances 
provide a general trend of their dependencies on A and and also indicate quantitatively 
the systematic uncertainties. 

We stress again, that HSD and UrQMD do not include explicitly a phase transition to the 
QGP. The expected enhanced fluctuations - attributed to the critical point and phase transition 
- can be observed experimentally on top of a monotonic and smooth 'hadronic background'. The 
most promising signature of the QCD critical point would be an observation of a non-monotonic 
dependence of the scaled variances on bombarding energy Ei a b for central A+A collisions with 
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fixed atomic mass number. In the fixed target SPS experiments the centrality selection in A+A 
collisions is defined by the number of the projectile participants. The measurements of cj_, uj + , 
and u c h are then preferable in the forward hemispheres. In this case the remaining small 
fluctuations of the number of target participants in the 1% most central collisions become 
even less important, as they contribute mainly to the particle fluctuations in the backward 
hemisphere. Our findings should be helpful for the optimal choice of collision systems and 
collision energies for the experimental search of the QCD critical point. 
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